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Abstract. We study a mean-field relativistic model which is able to describe 
both the behavior of finitely many spin- 1/2 particles like electrons and of the 
Dirac sea which is self-consistently polarized in the presence of the real par- 
ticles. The model is derived from the QED Hamiltonian in Coulomb gauge 
neglecting the photon field. All our results are non-perturbative and mathe- 
matically rigorous. 
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For heavy atoms, it is necessary to take relativistic effects into account. How- 
ever there is no equivalent of the well-known A'^-body (non-relativistic) Schrodinger 
theory involving the Dirac operator, because of its negative spectrum. The correct 
theory is Quantum Electrodynamics (QED). This theory has a remarkable predic- 
tive power but its description in terms of perturbation theory restricts its range 
of applicability. In fact a mathematically consistent formulation of the nonpertur- 
bative theory is still unknown. On the other hand, effective models deduced from 
nonrelativistic theories (like the Dirac-Hartree-Fock model [15]) suffer from incon- 
sistencies: for instance a ground state never minimizes the physical energy which 
is always unbounded from below. 

Here we study a variational model based on a physical energy which can be 
minimized to obtain the ground state in a chosen charge sector. Our model describes 
the behavior of a finite number of particles (electrons), coupled to that of the Dirac 
sea which can become polarized. Although it plays a minor role in the calculation of 
the Lamb-shift for the ordinary hydrogen atom (comparing to other electrodynamic 
phenomena), vacuum polarization is important for High-Z atoms j39| 146] and even 
plays a crucial role for muonic atoms [HI [21] . We show that the introduction of the 
vacuum in the model is the solution to deal with the negative energies of the Dirac 
operator and obtain a well-defined ground state. This was predicted by Chaix and 
Iracane in [51 page 3813]. 
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Our results are fully non-perturbative and mathematically rigorous. The cor- 
responding proofs are lengthy and therefore published elsewhere [231 IMl [Ml [21] ■ 
Unfortunately, we have not yet been able to include the photon field in the model 
for mathematical reasons, but a model with photons can be formally written fol- 
lowing our ideas. We emphasize that our goal is not to obtain all the QED effects 
accurately but rather to show how the introduction of the self-consistent vacuum 
changes dramatically the general properties of the model, leading to a well-defined 
variational theory. The fact that optimal states are found by a minimization prin- 
ciple is important for computational purposes and is essential for a justification of 
Relativistic Density Functional Theory (TSl [H] . 

Our methodology is as follows. We consider a Hartree-Fock type model in which 
particles interact through the Coulomb potential and with a kinetic energy given by 
the Dirac operator. Since we do not normal-order the underlying Hamiltonian, the 
kinetic energy is unbounded from below. However, we can as a first step construct 
the free Dirac sea by means of a thermodynamic limit. It is formally the minimizer 
of the Hartree-Fock energy. This state is not the usual sea of negative electrons 
of the free Dirac operator because all interactions between particles are taken into 
account, but it corresponds to filling negative energies of an effective mean- field 
translation-invariant operator. As a second step we introduce an external field 
potential and obtain a bounded-below energy by subtracting the (infinite) energy 
of the free self-consistent Dirac sea. In other words, we use the translation-invariant 
free vacuum as a reference and describe variations compared to it. We emphasize 
that this methodology is general and can be applied to other infinite quantum 
systems. It was used for the modelling of defects in crystals in [S]. 

1. Formal derivation of the model 

We start with the formal QED Hamiltonian written in Coulomb gauge, in the 
presence of an external electromagnetic potential (V, a), see |28[ [HI [13 [Ul [S] 

(1) H- ^ / *■(.) [„ ^ (-.V - Ai.) - a(x„ + „fl *(x, + / K(x).(x) dx 

, ff p{x)p{y) 



In this formula, ^'(x) is the second quantized field operator which annihilates an 
electron at x and satisfies the anticommutation relation 

The operator p(x) is the density operator defined by 

(3) ^^^^^^IKM^ 

(T=l 

where [a, b] — ab — ba. The operator Hf describes the kinetic energy of the photons: 

Hf = — I (|V X A(a;)p + \Et{x)\^) da; = - V / dk \k\a*Jk)ax{k) + Cte 
STra J ' 'a J^s 

(Cte indicates an infinite constant). The operators A{x) and Et{x) are the elec- 
tromagnetic field operators for the photons and a\{k) is the creation operator of a 
photon with momentum k and polarization A. 

In {V,a) is an external electromagnetic potential, for instance created by a 
set of nuclei. We use the notation 

L>° = -ia- V + m/3 



A MINIMIZATION METHOD FOR RELATIVISTIC ELECTRONS 



3 



for the Dirac operator. The constants m and a appearing in ^ are respectively the 
(bare) mass and (bare) Sommerfeld fine structure constant for the electron. The 
units are chosen such that h = c = 1. The Hamiltonian H^^° formally acts on the 
Fock space, 

where J-^ is the fermionic Fock space for the electrons and T-p\i is the bosonic Fock 
space for the photons. 

We emphasize that ([ij does not contain any normal-ordering or notion of (bare) 
electrons and positrons: ^(a;) can annihilate electrons of negative kinetic energy. 
The distinction between electrons and positrons should be a result of the theory 
and not an input. The commutator used in the formula ^ of p{x) is a kind of 
renormalization, independent of any reference. It is due to Heisenberg j28j (see 
also [m Eq. (96)]) and it is necessary for a covariant formulation of QED, see [47l 
Eq. (1.14)] and [131 Eq. (38)]. More precisely, the Hamiltonian H^-" possesses the 
interesting property of being invariant under charge conjugation since the following 
relations hold formally 

where is the charge conjugation operator acting on the Fock space. 

In our study of the QED Hamiltonian H^'", we shall make two approximations: 

• we neglect photons and assume there is no external magnetic field, a = 0; 

• we work in a mean-field theory, i.e. we restrict the Hamiltonian to Hartree- 
Fock states. 

These approximations are of a different importance. Neglecting photons is of course 
a very rough approximation as it will forbid us to describe important physical effects 
occurring in QED like the self-energies of the electrons, the biggest contribution to 
the Lamb shift. But we do that only for mathematical reasons: we were not yet able 
to extend most of the results presented below when photons are taken into account. 
Formally^ a large part of our study is exactly the same with photons (when they 
are treated by a mean- field procedure) . We hope to come back to this point in the 
near future. 

The second approximation which we make by restricting ourselves to Hartree- 
Fock states is more fundamental and many of our results are specific to this case. 
Nevertheless, some of our general ideas may be applicable to the full QED model. 

Let us recall that the electronic one-body density matrix (two point function) of 
any electronic state \ E J-c is defined as 

Pix,y),,^, = {n\-^*{x),-^iy),,\n). 

In view of ([3]), it is natural to introduce a renormalized one-body density matrix 



7(2;, y)^,^' = 

By ©, we obtain the simple relation 



' 2 

where / is the identity operator. Electronic Hartree-Fock states form a subset 
{|J7p)} C J-Q of states which are completely determined by their density matrix P 
(or equivalently by their renormalized density matrix 7 = P — 1/2). Recall that if 
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is a Hartree-Fock states with N occupied orbitals (pi, ...,(pN, then the associated 
density matrix P is just the orthogonal projector on Span(93i, ...,ipn): 

N 
1=1 

For a formal Hartree-Fock state with infinitely many occupied orbitals 
we also obtain 



Hence 



„ / P- P 

7 = P-7T = 



\ occ unocc / 



2 2 

The associated density of charge is formally given by 

(4) p,{x) = {n\p{x)\n) = ^ ( E - E \^^(^)n ■ 

\ occ unocc / 

Now we can compute the energy of any state \ilp) (Xi |0) where \ ftp) is a Hartree- 
Fock state in J^c and |0) S J^ph is the photonic vacuum. We obtain 

(o| (g) {np\m^'°\np) ® |o) = s^^iP - 1/2) + cte 

where Cte is an infinite constant and 



(5) £:^f(7) - triD'l) + / V{x)p^{x) dx 



a 



p^{x)p^{y) a // ItOm/)^ , , 
■dxdy—— // — —dxdy. 



2JJ \x-y\ 2JJ \x-y\ 

The reader can recognize the well-known Hartree-Fock energy, but applied to the 
renormalized density matrix 7 = P — 1/2 instead of the usual density matrix P. 
The last two terms of the first line are respectively the kinetic energy and the 
interaction energy of the electrons with the external potential V. In the second line 
appear respectively the so-called direct and exchange terms. In Relativistic Density 
Functional Theory [151 116j , the exchange term is approximated by a function of p^ 
and its derivatives only. 

Any stationary point of the above energy satisfies the first order equation (written 
in terms of the usual density matrix P — j + 1/2) 

[P, Fp_i/2] - 
where -Fp_//2 is the Fock operator 

^''-'/--^° + '^ + °''i'--'/'i'R-° '^J-T''' ' 

For a minimizer (in a chosen charge sector) , one will have the more precise equation 

P = X(-oo,m] {Pp-1/2) 

where /it is a Fermi level and X(-oo,/i] (^) is a mathematical notation for the spectral 
projector of A corresponding to filling all energies < /z. Saying differently, one 
obtains a Hartree-Fock state with infinitely many occupied orbitals, all having an 
energy < p. We shall give a precise interpretation of this equation later on. 

It is time to worry about the mathematical meaning of the formulas we have 
formally derived up to now, in particular the definition of the energy ([5]). Unfortu- 
nately, the latter does not make any sense for the following reason: when P is an 
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orthogonal projector (as this is the case for HF states), 7 = P — I /2 is never a com- 
pact operator in an infinite dimension space. Hence none of the terms appearing 
in ((5]) has a clear mathematical meaning. Formally, one has £^p{P — 1/2) = —00 
for any density matrix P. 

In [26j . we proposed to overcome this difficulty in the following way: we restrict 
the whole system to a box of size L with periodic boundary conditions and an 
ultraviolet cut-off A in the Fourier domain. Then all the above formulas make 
perfectly sense because we are in a finite-dimensional setting. In particular one can 
define minimizers of the HF energy with or without the external field V, with or 
without a charge constraint. Then, we look at the limit of the minimizer in the 
considered class when the size of the box grows, L ^ 00, but the cut-off A stays 
fixed. The limit (if it exists) is the formal minimizer of the unbounded below energy 

C-jjp. 

Actually we shall essentially use this method to define the free vacuum (the 
global minimizer of when V = 0). Once the free vacuum has been found, we 
formally subtract its (infinite) energy to the expression (O and obtain a well-defined 
bounded below energy called Bogoliubov-Dirac-Fock and which is related to a work 
of Chaix and Iracane [1] . 

Notice the ultraviolet cut-off A is fixed during the whole study. It is only at the 
very end that we can tackle the difficult task to remove it by renormalization. We 
shall also discuss the appearence of the Landau pole. 

We explain all that in details in the next sections. 

2. Restriction of the system to a box and definition of the free 

VACUUM 

Let us consider a box of size L, Cl '■= \—L/2;L/2Y and limit the system to 
this box, with periodic boundary conditions. For simplicity, we also periodize the 
Coulomb potential and introduce 

(6) ^-(-) = z^ E 1^^"^+^^^ 

Vfee^\{0} 

where c is chosen such that G > 0. Furthermore, we add an ultraviolet cut-off A, 
i.e. we choose as one-body space the finite-dimensional 



f)i := Span ( e**--^ 



fce-^, |fc|<A 



The periodic Hartree-Fock energy (without external field V) is then defined as 



(7) £:°(7l) = triD'^lL) + | / / P^dx)GL{x - y)p-,Mdxdy 



\lL{x,y)\ Gl{x -y)dxdy. 



This expression is well defined for all renormalized density matrices 7l acting on 
the one-body space and satisfying the condition that "fL+ 1/2 is an orthogonal 
projector. Indeed, following a method of Lieb [37], we can even relax this condition 
and work under the assumption that 

(8) - 1/2 <-iL< 1/2. 

It is possible to define the QED Hamiltonian without photons in the box in the 
same way, see [26j. Notice the fermionic Fock space built on the one-body space 
is also finite-dimensional. 
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The minimizing problem defining the free HF vacuum in the box reads 

El := inf SUjl). 

-//2<7i<//2 

It was shown in [26l Thm 2.7] that for L ^ 1 and < a < 4/7r this problem admits 
a unique minimizer 7°, which has several interesting properties. First it takes the 
form jI = Pj^ — //2 where P^^ is an orthogonal projector acting on Sjj(, hence the 
relaxation ^ does not change the minimum. Then, 7^ is a translation-invariant 
operator, meaning that it is a multiplication operator in the Fourier domain, 7^ = 
7° (fc). It can also be proved that the associated density of charge vanishes, p^i^ = 0. 
This comes from the fact that 7° has a very special form which we do not detail as 
we are more interested in the properties of the limit of 7^ as L ~^ 00. 
Indeed, it was shown in Thm 2.7] that 

uniformly as functions of the Fourier variable and that 

as L — > 00. The operator 7*^ is the density matrix of the free Hartree-Fock vacuum 
in the whole space (with the ultraviolet cut-off A), which formally minimizes the 
no- photon QED Hartree-Fock energy £^p in spite of the fact that its energy is —00. 

We now describe the interesting properties of 7*^, which were proved in [26|, Thm 
2.2] and [35]. First 7° = I'^ip) is a translation-invariant operator acting on the 
one-body space 



{/gL2(R3,c4), Supp(/) Ci3(0,A)} 



of functions whose Fourier transform is supported in the ball of radius A (the natural 
"limit" of SjjO- One has 7" = — 1/2 where is an orthogonal projector which 
satisfies the following SCF equation: 

=X(-oo;0](2?°), 



(10) S -no _ no „,{vl-i/2)ix-y) 



\x-y\ 

The operator I?" is the self-consistent (SCF) Fock operator of the free vacuum V'^ 
It was shown that it takes the following special form 

2?°(P) = .9i(bl)a-p + .9o(bl)/3 

with 

l<.i(bl)<^, 
m 

hence 



(11) \D'{p)\ < \v"ip)\ = Vg^mW+M?, 

i.e. the gap of T)'^ is bigger than the one of the original Dirac operator Notice 
(fTU]) corresponds to the usual Dirac's picture that the free vacuum is a Hartree- 
Fock state occupying all the negative energies of a Dirac-type operator. If a = 
(no interaction), then we get the original picture = P° := X(-oo;0] (^°)i but in 
general V°_ ^ P^. 



-po 

Figure 1. The free vacuum fills the negative energies of the 
SCF Fock operator 
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Notice ([TUI can be rewritten in terms of 7° in the form 
(12) 

2|I?"(p)| 2^g^{\p\YW+go{\p\Y 2./^:M?W^M[f 

In QED, the Feynman propagator at equal times 

Sf{x, y; = ty) i7(x, y)l3 

is often expressed using the Kallen-Lehmann representation [3T1 E] , based on 
relativistic invariances. Although our model is not fully relativistically invariant (we 
discard photons and use an ultraviolet cut-off A) and is only defined in the mean- 
field approximation, our solution (|12p has exactly the form which may be derived 
from the Kallen-Lehmann representation for the equal time propagator. In four- 
dimensional full QED, a self-consistent equation similar to pop is well-known and 
used. These so-called Schwinger-Dyson equations [48l [M] have been approximately 
solved for the free vacuum case first by Landau et al. in (TJ [2] , and then by many 
authors (see, e.g., [301 1211 13] ) ■ Equation pH)) has already been studied by Lieb and 
Siedentop in [3S] in a different setting. 
We notice that 

p^o = 0. 

This is indeed a consequence of Formula (fT2|l : one has Cj^C^^ = —7° where C 
is the charge conjugation operator. Hence any negative energy state of 7*^ can be 
associated to a positive energy state obtained by charge conjugation. The result 
follows from In mathematical terms, p^o{x) — tic4,j'^(x,x) = 0, the Dirac 
matrices being trace-less. 

We want to mention a last interesting property of 7*^: it is indeed the unique 
minimizer of the energy per unit volume defined by 
(13) 

r(7) = ^ / trc.[i.°(p)7(.)]rfp- ^ // 

{^■^r J\p\<A i^T^r J J\p\^q\<A \p-qr 



where we recall that A is the ultraviolet cut-off. This property can be used for the 
numerical computation of the free vacuum 7^^. Lastly, we have that the energy per 
unit volume of the free vacuum is 

e = inf r(7), 

-I/2<1<I/2 

the limit appearing in ^ as proved in j26i Thm 2.7]. 



3. BOGOLIUBOV-DlRAC-FOCK THEORY 

If we summarize, using a thermodynamic limit we have been able to define the 
free vacuum which is the unique minimizer of fyp, and despite the fact that its 
energy is —00. The free vacuum is the negative Dirac sea of an SCF translation- 
invariant Dirac Fock operator D^. Now we use this vacuum as a reference and 
subtract its (infinite) energy to the original HE energy, in order to obtain a bounded- 
below function. Formally, this gives for any state 7 = P — 1/2 the so-called 
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Bogoliubov-Dirac-Fock (BDF) energy 110] 

^bdf(^ - 'P-) - (0| ® (rip|H^'°|f)p) ® |0) - (r!o|H"'°|r!o) 
= £^p(P-//2)-£°p(pO -//2) 

= trX>"(P-7'°)+ / V{x)pp^Tyo{x)dx 



- 1 // "'t-r - ? // 

where |r2o) = \^-po ) (g) |0) is the no-photon HF free vacuum in Fock space found in 
the previous section. In (|14p we have used that p-po_jy2 = and recognized the 
formula ([13 of X>°. 

The BDF energy measures the energy of any state 7 compared to the (infinite) 
energy of the free vacuum 7°. Also Q = P — describes the variations counted 
with respect to the free Dirac sea. The BDF energy was first introduced by Chaix 
and Iracane [9] but with 7-"° and V'^ replaced by P° and D°. It was first mathe- 
matically studied in [S] . Chaix and Iracane obtained their energy by imposing from 
the beginning a normal-ordering on the QED Hamiltonian, taking as definition of 
positrons and electrons the ones given by the decomposition induced by P^. If 
a = our model is equivalent to the one of Chaix-Iracane, but it is not when 
a 7^ 0. It seems that normal ordering is only fully relevant for the description of 
non interacting systems. 

Once again the above formal computation (jl4p can be justified by a thermody- 
namic limit. We will show that the last expression of S^j^p is well-defined math- 
ematically and we will be able to find minimizers of this energy. We can prove 
that any sequence of minimizers in boxes will converge to these states in the ther- 
modynamic limit L —t 00 but we do not give more details and refer to |261 Thm 
2.9]. 

We now explain how it is possible to give a mathematical meaning to the last 
expression of (|14p . Some details which may appear as mathematical technicali- 
ties will later reveal to be crucial for renormalization, hence related to important 
physical properties. We recall that an operator Q is said to be trace-class when 
Till {'^iWQ*QWi) < 00 in some orthonormal basis {fi) of the one-body space. 
Then Xr{Q) = Ti {'■Pi\QWi) is well-defined and does not depend on the chosen ba- 
sis. In principle it is possible that the series {'^i\Q\^i) converges for one specific 
basis even if the operator is not trace-class. This will be the case for our operator 

Given an operator Q, we define :— V^QV^, where e, e' G {±} and := 
1 — V?.. We say that an operator is P^ -trace class if Q"*"^ and Q are trace-class 
and we define 

tro(Q) tr(g++) + tr(Q-) - ^ {vt\QWt) + J2 i^^lQl^i) 

i i 

for any chosen basis {ff ) U {ifiY) adapted to the decomposition induced by . Of 
course if Q is trace-class then it is also P^-trace class but the converse is not true. 

Now we remark that when tr(P — V^)'^ < 00 for a projector P, the operator 
Q = P — is automatically P° -trace class. The reason is that 

(P - P° )2 = (P - P° )++ - (P - P° )~. 

Additionally tro(P — P" ) is always an integer as proved in [531 Lemma 2] and [1]. 
We interpret etro(P — V^) as the charge of the state P (measured with respect to 
the free vacuum). Notice the condition tr(P — V?_Y < 00 is a classical requirement 
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of the Shale-Stinespring Theorem fST which guarantees equivalence of Fock space 
representations . 

Now we notice that when P — is V^l-tiace class, 

tro(I?°(P - Vl)) = tr(|I?0|((P - Vl)++ - (P - Vl) — )) = tr \V"\{P - Vl)'^ > 0, 

i.e. the kinetic energy is non negative and well defined when tr(P — Vl)^ < oo. 
Using Kato's inequality \x\^^ < (7r/2)|p| and (fTTjl we infer, following [5], 

(15) f // '^^'iJl^^^'^^'^ x^. < ^tr(H(P-P°)^) < ^tro(P"(P-P^)), 

i.e. the last term of (fT4|) is also well-defined. 
Now we assume that 

V = —av * 1 — : 

is the electrostatic potential created by a set of extended nuclei with (fastly decaying 
and smooth) total density j/, J v = Z . We define the BDF energy of Q = P — P'^ 
by 

£lMQ) ■■= ir,{V^Q)-aD{pQ,v) + |i?(pQ,PQ) - III M^^dxdy 
where 

is the so-called Coulomb scalar product. 

It was proved in |25l Lemma 1] that when trQ^ < oo and Q is P° -trace 
class, then pq is a well-defined function which is squared-integrable and satisfies 
D{pQt Pq) < oo, hence ^^j^^iQ) well-defined by ([TS]). Additionally we have when 
< a < 4/7r by p3)) and using that D{-,-) is a scalar product 

^BDF(^-^-)>-f^(^,i^)>-00, 

hence the BDF energy is bounded from below. 

After these mathematical details, we are now able to minimize the BDF energy. 
We can either look for a global minimizer which will be interpreted as the polarized 
vacuum in the presence of the external potential V, or for a minimizer with a charge 
constraint 

tro{P-Vl)=N 

which will usually represent the state of N electrons coupled to the self-consistent 
polarized vacuum. We detail the two situations in the next sections. In both cases, 
the obtained minimizer will be P° -trace class but not trace class (except when 
V = 0), which will be related to renormalization as we will explain later. 

4. The Polarized Vacuum 

The polarized vacuum is by definition the state of lowest QED energy in the 
Hartree-Fock no-photon class. By (HU, it is also the state of lowest BDF energy. 
Hence we consider the following minimization problem 

E"" :=inf£^Dp(P-pO) 

where the minimization is done over all orthogonal projectors P acting on Sja such 
that P — P° is P° -trace class. As before the constraint on P can be relaxed 
following Lieb [37] and replaced by the convex constraint < P < /, but we do not 
detail this here. 
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It was proved in [231 I24j that a minimizer Pvac exists and that it solves the 
self-consistent equation 

J ^Vac = X(-oo;0] (-P'Pvac-//2) J 

Hence one more time the vacuum Pvac corresponds to filling negative energies of a 
self-consistent Fock operator. Notice 

(17) Pp__//2 = Po + 0(a2) Fo = DO + V. 

^ ^ -^o — o—or A — — I 



vac 

Figure 2. The polarized vacuum Pvac in the presence of V fills 
the negative energies of the SCF Fock operator Pp^^^_//2- 

In general, one could have to create electron-positron pairs if one wants to deform 
V'^ into the polarized vacuum Pvac- But when V is not too strong it was proved in 
that Pvac is unique and neutral: 

tro(Pvac-P-) = 0. 

In this case the vacuum Pvac only contains virtual electron-positron pairs compared 
to P° , see [221 Appendix] . 

In the right hand side of (|17p . P[p^^^-i/2]j represents the vacuum polarization 
density, which is self-consistently created by the external potential V. Notice that 
one of the highlights our procedure is that although the reference V'^ appears in 
the functional f , the equation (fT6|) is independent of P*^ , showing that the free 
vacuum energy serves just as helpful device. 



5. Atoms and molecules 

For the study of common physical systems like atoms or molecules we have 
to consider the minimization of the BDF energy in charge sectors, that is to say 
imposing a constraint of the type 

tro(P-P*l) =iV 

where iV G Z. Of course we cannot impose the number of particles but if V is not 
too strong and TV > 0, this will provide a system of N electrons interacting with 
the vacuum. Hence we introduce the following minimization problem 

(18) E''{N):= inf E^MP-rl) 

where as before P is assumed to be an orthogonal projector such that P — P'l is Pi- 
trace class. It is not expected that a minimizer will always exist. If for instance N is 
too large compared to the number Z of nuclei, the system will certainly be unstable. 
On the other hand if Z is too large, pairs could be created, which complicates the 
description of the system. In 25], it was proved that when the following binding 
conditions hold true 

(19) E^{N)<E^(N~k) + EO{k) VfceZ\{0}, 

then a minimizer exists for E^ (N). The binding condition (|19p was proved to hold 
true in when for instance < N < Z and a <C 1 (non relativistic limit). 
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Here we assume that there is a minimizer P and that V is not too strong. Then 
it was proved in [25] that P solves the SCF equation 

, . } P = X(-oo;m] {Fp-1/2) , 

\ Fp-i/2 ^D^ + a(p[P-//2] - * R - \Ly\ 

where /i is a Fermi level (a Lagrange multiplier due to the charge constraint). We 
can write 

P Pva.c ^" Pel 

where 

N 

^'vac = X(-oo;0] {Fp-i/2) and Pol = X(0;m1 {.Fp-1/2) = X! 

i=l 

with 

for all eigenvalues A,; < The orbitals (vi)fLi describe the Hartree-Fock state of 
the N electrons whereas Pvac describes the SCF polarized vacuum in the presence 
of the external field V and the N electrons. 



^1— ocH^ zn 



1^ 



Figure 3. Decomposition of the system 'vacuum + TV electrons' 
for the solution P in the iVth charge sector. 

We notice that the decomposition of the state P into N electrons and the po- 
larized vacuum can be made unambiguous because P satisfies the SCF equation 
(PO)) . For a general state P satisfying tro(P — V^) — N, there is no canonical 
decomposition between real and virtual particles. 

Now we remark that 

(21) Pp-//2 = ^ci + aP[p,,,-//2] * f-r - " 1 _ I = Fci + 0{a ) 

fI p y\ 

where 

Pol = P + a(pp„i - ^) * - r 

\x\ \x-y\ 

is the usual Dirac-Fock operator for N relativistic electrons. Hence we deduce that 
the (fii 's solve the usual Dirac-Fock equations [IHl [TT] , perturbed by the SCF vacuum 
polarization potentials. An essential feature is of course that these equations have 
been obtained by a minimization procedure, contrarily to the usual Dirac-Fock case. 

Notice the Dirac-Fock model is not obtained as a variational approximation of 
the BDF model. But the Dirac-Fock equations are an approximation of the BDF 
equations. This was first noted by Chaix and Iracane in ^9J. 

6. Time-dependent equation 

The time-dependent equation corresponding to our model could also be useful, 
in particular for the study of spontaneous pair creation which is usually formulated 
as an adiabatic theory on the evolution equation [JDl |33] . It reads 

tPit) = [FlP^t)-I/2],Pit)] 
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where we choose as initial condition an orthogonal projector P{0) such that tr(P(0)— 
7-"° )^ < oo. It was proved in 27J that this equation admits a global-in-time solution 
P{t), i G M, which has a constant BDF energy and charge: 

7. Renormalization 

In regular QED, the divergences of the (appropriately defined) physical measur- 
able quantities are usually eliminated by means of a mass and a charge renormal- 
ization. The main idea is to assume that the parameters a and m appearing in 
the theory are indeed bare parameters which are not physically observable. The 
physical parameters are assumed to be functions of a, m and the cut-off A 

"ph = aph(a,TO, A), TOph = mph(a,m, A) 

and equal the physical values obtained in experiment. These functions should be 
inverted in order to express the unknown bare quantities in term of the physical 
quantities 

(22) a = Q;(aph,mph, A), m = m(aph, mph, A). 

Using these functions, one expects to remove (in some sense that needs to be pre- 
cised) all divergences from physically measurable quantities. 

Mass and charge renormalization however does not remove all divergences in the 
theory. Certain quantities, e.g. the bare Feynman propagator Sp (either at equal 
times or at general space time points), are still divergent. The expectation is that 
all these divergences cancel in physically measurable quantities and that they are 
therefore of no real relevance in formulating the theory. 

Although there is no real need to do this, it is often convenient to introduce 
a renormalization of the bare Feynman propagator Sp- This is referred to as a 
wavefunction renormalization. In full QED [14j it is claimed that the divergence 
in the Feynman propagator may be removed by a multiplicative renormalization 
and that the renormalized propagator has the same pole near mass shell in 4- 
momentum space as a free propagator corresponding to a particle with the correct 
physical mass. 

Note that in practice, this theoretical renormalization procedure is always used 
to justifying the dropping of the divergent terms obtained at each order of the 
perturbation theory _14J. For this fact to be true, it is particularly important that 
renormalization can be expressed by means of multiplicative parameters in front of 
the different propagators JT^. 

In Hartree-Fock QED, it is not clear at all if the usual renormalization program of 
QED can be applied, especially when photons are not included. In [44, p. 194-195], 
it is argued that mass and charge renormalization is alone not enough to completely 
remove the divergences of the HF theory by means of multiplicative parameters. 

In any case, the physical mass and charge have to be identified within the model. 
We propose the following definitions. The physical mass is just the lowest energy 
of a free electron, hence 

(23) mph(a,m,A) :=£;"(1) 

which was defined in (fTH]) . 

To define the physical coupling constant, we consider an extended nucleus of 
density v, J v — Z , and put it in the vacuum. Let Qvac = -Pvac ~ P- be the 
polarized vacuum solution of (|16p . We assume that v is not too strong such that 
the vacuum stays neutral, tro((3vac) =0. Of course in reality it is impossible to 
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distinguish the nucleus from the vacuum and the charge which is observed far way 
from the nucleus is just 



e ^ - / PQvao 

(provided pq^.^^ is an function). Hence we may define 

(24) aphys(a, m, A) := a ^1 - Z^^ J pg^^^^ . 

If the above formula still depends on Z, one can take the limit as Z — s- 0. 

It is very important to realize that charge renormalization is based on the fact 
that the operator Qvac is not trace-class. If it were trace-class, one would of course 
have tro((5vac) — = / pq^^^, hence aph — a. Therefore, the mathematical diffi- 
culty that a minimizer of the BDF energy is never trace-class (except when j/ = 0) 
is the origin of charge renormalization. Also this shows that in a finite dimensional 
space (for computational purpose for instance), renormalization is certainly more 
involved as all operators are trace class. 

Both (|23|) and (|24|) would define TOph and ctph as extremely complicated non- 
linear functions of a, m and A. A challenging task is to study the finiteness of 
measurable quantities like for instance the energy of an electron in the presence of 
an external field i?^(l), when oiph and mph are fixed to be the observed physical 
quantities. We do not know if this is possible when photons are not taken into 
account. 

It is however possible to completely solve the above program for a (further) 
simplified model called the reduced Hartree-Fock, as was done in [23j. We explain 
that now. 

The reduced HP model is just obtained by neglecting the exchange term in the 
HF energy Q 

(25) CF(7) = tr(i?%) + j V{x)p,{x)dx + ^ J J Eli^kM da: dy . 

This is natural as the exchange term is usually treated together with a term describ- 
ing the interaction with the photon field to form the standard electron self-energy 
that is a subject of the mass renormalization. 

The so-obtained model is much simpler than the HP model as the energy is now 
a convex function of 7. All what we have said concerning the case with exchange 
term can be extended to this simplified model. The free vacuum is even a simpler 
object as in Pquation (jTU]) only the exchange term created a self-consistent field. 
Hence we obtain 

Vl=P° and I?" = i?". 
The reduced Bogoliubov-Dirac-Fock (rBDF) energy then reads [24] 

(26) Cdf(Q) - tro(i^"Q) - aD{pQ,v) + |d(pq,pq). 

It can easily be shown that for a free electron in the vacuum 25, Lemma 3] 

inf £,%p(P-P^) = m, 

tro(P-P°) = l 

i.e. TTiph = m and there is no mass renormalization for the reduced BDF model. 

Consider now a small external density v, ^ v — Z and let Qvac be the associated 
polarized vacuum, with density pvac PQ^^a- The SCP equation satisfied by Qvac 
reads 

(27) Qvac-X(-oo;0](J^)-^- 
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where 



F = D° + a{p^ 



We expand ^ in powers of a, using that ^ o'(-F') when v is small enough. 
We can use the resolvent representation [32l Section VI, Lemma 5.6] to derive the 
self-consistent equation for the density pvac 



1 

(28) Pvac{x)^-— dr/Trc 



DO + a{pQ - i^) * |i| + D0+ iij 



{x,x). 



Applying the resolvent equation 
1 1 
A 



1 1 



-2 -B-B- 



, 1 „ 1 „ 1 „ 1 

a^-B-B-B- 

A A A A-aB 



A~aB A A A AAA 

and using Furry 's Theorem |19j . telling us that the corresponding a^-term with two 
potentials vanish, we obtain 

(29) pvac = a.Fi [pvac - H + ^3 [^(Pvac ~ v)] , 



ZTT 



1 1 
— p * 



x\ D° + ii'/ * |a;| DO + irj"^ \x\ D^ 



p * 



{x,x). 



As realized first by Dirac [II1II2] and Heisenberg [28j, cf. also [20^, the term Fi [p] 
plays a particular role since it is logarithmically ultraviolet divergent. Following, 
e.g., Pauli-Rose [H], one evaluates in Fourier representation 

F,[p]{k)^-p{k)BAik), 

where [42l Eq. (5)-(9)] Ba(A:) =Ba- Ca(A:), with 
(30) 

Ba = Ba(0) = - 

TT 

and 



/T+a2 Z' 



zV3 



1-z^ 



:^log(A)-^ + :^log2 + 0(l/A2). 
37r 97r 37r 



(31) 



lim CA(fc) = C(fc) 

A^oo 



1 

2^ 



dx{l - x^) log[l + k^{l - a;2)/4]. 



which was first calculated by Serber and Uehling [451 EQ] . 

We can now compute the physical coupling constant. First we rewrite ((^^ in 
Fourier space as 



(32) (1 + aSA)pvac(fc) = aBfJ){k) + aCA{k){p,,, - P)(fc) + FsHp, 
Assuming that pvac G L^(M^) and taking fc = 0, we find 

aBj^Z 



iy)]{k). 



^0 



l + aBA 

where we have used that Ca(0) = -F3[a(pvac — t^)](0) = 0. Hence by ([24]) we find 
(33) 



1 + a^A 



It follows that necessarily aph-BA < 1- We emphasize that although in the literature 
the expression of ctph is sometimes expanded to get aph — a{l — aB^) leading to 
the condition aB\ < 1, the real constraint indeed applies to the physically observed 
Q!ph and not the bare one. 
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We now show how to renormaUze the SCF equation using (|33p . Denote p 
Pvac — I' the total (observable) density, then ()29p can be rewritten in terms of p 



(34) 
and 
(35) 



ap 



a^BAP + a^CA{k)p + aFsiap] 



ap 



-CA{k)ap 



a 



-Faiap]. 



1 + uBa 1 + aBA 1 + aBA 

To perform our renormalization scheme we fix as physical (renormalized) objects 
ttphPph = ctp. Notice the renormalization of the density p is similar to a wavefunc- 
tion renormalization of the (equal time) Feynman propagator as explained above. 
We can rewrite the self-consistent equation as 



(36) 



^phPph 



-aphi^ + alhCA{k)pph 



aph-F3[Q;phPph], 



independently of the bare a. Notice that equation (I36p satisfied by Ofph/Oph is exactly 
the same as equation (j34[) satisfied by ap, but with the logarithmically divergent 
term a^BAP dropped. Therefore, as usual in QED the charge renormalization 
allows to simply justify the dropping of the divergent terms in the self-consistent 
equation. In practice [3^, one would solve with ofph — 1/137 and with Ca(A:) 
replaced by its hmit C{k). 

Returning to the effective Hamiltonian F = £'*' + Q;(/3vac^'^)*l/|a;| and inserting 
pop , i.e. expressing in terms of the physical objects, we obtain 

(37) D" + aphpph * ^ D° - aphi^ * + Vcs , 



with 



CA{k)pp\y{k) + aph^3(aph/3ph) 



fc2 



the effective self-consistent potential, where denotes the inverse Fourier trans- 
form. Notice, this equation is valid for any strength of the external potential. 
However, expanding pph in aph, we obtain to lowest order in aph 



off 



2 ^ 

"ph 3 



a 



Ph 

37r 



CA{k)Hk) 
k^ 

1)1/2 



(x) 



2 

<2 



1 



Kx') 



the Uehling potential [7]. 

8. The Landau Pole 
We notice that can be written as 

Oph 

a = . 

I - ttphi^A 

The fact that the denominator can go to zero is usually called the Landau pole. 
Also we see that 

(38) aph^A < 1 

which proves that aph when A oo, independently of a. 

In [24,, Thm 2] , it was proved that for a fixed (and not too strong) external field 
V = —av * jiy, the unique polarized vacuum Pa of the reduced BDF model satisfies 

lim tr (Pa — P")^ ~ and lim D(pp po — v,pp pa — v) — 0. 

In words, when A — > oo, the vacuum polarization density totally cancels the 
external density for pp^^po v. But since Pa — P- ^ 0, this means that in the 
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limit A oo, Pa — P- and its associated density become independent. Therefore, 
the minimization without cut-off makes no sense both from a mathematical and 
physical point of view. Indeed all this easily implies that when no cut-off is imposed 
and when v ^ 0, the infimum of the reduced BDF functional is not attained. In 
physics, this "nullification" of the theory as the cut-off A diverges has been first 
suggested by Landau et al. [Ml [331 HI 133 and later studied by Pomeranchuk et al. 

m- 

We notice that with the usual value aph — j^y, ([35]) leads to the physical bound 
A < 10^*° (in units of mc^). 

9. Conclusion 

We have presented a model which is obtained as the mean-field approximation of 
no-photon QED. We believe that the Hartree-Fock approximation is an interesting 
model as it possesses already many peculiarities of the full QED and it is much 
simpler to handle. In particular, optimized states always correspond to filling the 
spectrum of a one-body operator up to some Fermi level /i, which corresponds to 
the original interpretation of Dirac. 

The main advantage of this model is that it is variational: states can be found 
by minimizing an energy, contrarily to the usual relativistic effective models used 
for instance in Quantum Chemistry. This provides a better interpretation of the 
optimal states. Also the model provides a justification of the Dirac-Fock equa- 
tions, which are seen as a 0{a^) approximation of a set of equations obtained by 
minimization. 

Another advantage of the model is that it is nonperturbative: the only constraint 
to have a globally stable model is that < a < 4/7r. The equations are quite simple 
and renormalization can be done non perturbatively to all orders (at least when 
the exchange term is neglected). 

The main idea in the derivation of our model was first to define the SCF free 
vacuum by a thermodynamic limit, and then to subtract its infinite energy in 
order to get a bounded-below function. This method replaces the usual normal- 
ordering which can only be used for non interacting systems. In principle the same 
method could be used for the full QED. But probably it is not possible to express 
the difference between the energy of the considered state and the one of the free 
vacuum in a simple way. 

We have neglected photons but in principle one could take the photon field into 
account. The mathematical study of such a theory remains to be done. 

Acknowledgment. M.L. and E.S. acknowledge support from the ANR project 
"ACCQUAREL" of the French ministry of research. 
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